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We consider monojet searches at the Large Hadron Collider (LHC) for spin-1 dark matter that 
interacts with quarks through a contact operator. If the dark matter particles are produced with 
longitudinal polarizations, then the production matrix element is enhanced by factors of the energy. 
We show that this particularly effective search strategy can test models for which the energy sup¬ 
pression scale of the operator is as large as 10® TeV. As such, these searches can probe a large class 
of models for which the contact-operator approximation is valid. We find that for contact operators 
that permit velocity-independent dark matter-nucleon scattering, LHC monojet searches for spin-1 
dark matter are competitive with or far surpass direct-detection searches depending on whether the 
scattering is spin independent or spin dependent, respectively. 


I. INTRODUCTION 

Broadly speaking, there are two search strategies for dark matter (DM) at hadron colliders: cascade decay searches 
and direct-production searches. The goal of cascade decay searches is to produce a heavy QCD-charged particle whose 
decays are required by symmetry to produce dark matter, and the prototypical examples are supersymmetry searches 
for squark or gluino production (see, for example, m)- The goal of direct production searches is to produce dark 
matter through a four-point (not necessarily contact) interaction that couples two dark matter particles directly to 
partons. The prototypical examples of this strategy are “mono-anything” searches [iUIO] in which one searches for a 
single jet, photon, W, or Z radiated from the partons, against which the DM particles can recoil. 

Cascade decay searches are most sensitive to the properties of the heavy particles that are initially produced, 
whereas it is typically more difficult to indirectly determine the properties of the DM decay product. On the other 
hand, although direct searches may have the disadvantage of not utilizing a strongly coupled production process, they 
can more directly probe the properties of the DM particle and its interactions with the Standard Model (SM). Direct 
searches have been used extensively at the LHC for the case where dark matter is a spin-0 or spin-1/2 particle [TTVllhj . 
The goal of this work is to discuss the unique signatures arising from the case of spin-1 dark matter (which we denote 

by 5)- _ 

A characteristic feature of vector dark matter is that the longitudinal polarization vector scales as ~ E/m, implying 
that the production matrix element can receive a significant enhancement in the region of phase space where the dark 
matter is relativistic and either one or both particles are longitudinally polarized m- The large enhancement of the 
matrix element results in large cross sections for the process pp —>■ B -I-jet, yielding greater sensitivity to the dark 
matter-quark coupling. 

The plan of this paper is as follows. In Section II, we tabulate the set of contact operators that we use to model 
dark matter-quark interactions, and determine the energy scaling of the DM production cross section that would be 
expected from considerations of angular momentum, charge conjugation, and parity. In Section III we discuss the 
limits on the applicability of the contact-operator approximation that arise from unitarity of the scattering amplitude. 
In Section IV we present the exclusion bounds on these models imposed by results of the 8 TeV LHC run, the future 
sensitivity that may be expected from the 14 TeV LHC run, and the implications for dark matter-nucleon scattering. 
We conclude with a discussion of our results in Section V. 


II. OPERATORS AND MATRIX ELEMENTS 

For simplicity we work under the assumption that spin-1 dark matter interacts with quarks through effective contact 
operators. The most general set of contact operators that we can write, up to dimension 6, is given in Table |Tj Note 
that we have assumed that H is a complex vector field. If H is a real field {B = 5^), then operators V(3 — 6) and 
V(7 — 10)_ vanish identically. The remaining operators must be scaled by a factor of 1/2 if B is real. We follow 
notation similar to [m, wherein the operators V(7 - 10) are herein referred to as V(7 - 10)+ Q 


^ Note, Ref. [17| did not consider the operators V(7— 10) — . For completeness, we redo much of the analysis of Ref. [17| for these operators 
in Appendix C. 
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Operator 

Structure 

Dim D 

VI 

{llA)BlB^qq 

5 

V2 

{l/K)lBlB'^q'i^q 

5 

V3 

(l/2A2)*(Bta^B- _ B''d^Bl)q-i^q 

6 

V4 

{l/2A^)l{Bld^B'' - B-'d^Bt)qj>^^^q 

6 

V5 

{l/A)lBlB,qa>^-'q 

5 

V6 


5 

V7+ 

{l/2A^){Btd-'B^ + Bod-'Bl)q-f>^q 

6 

V7- 

{l/2AyiBtd''B^ - B.d-'Bl)q-f>^q 

6 

V8+ 

il/2A^){Bld-'B^ + B,d''Bl)qy>^y^q 

6 

V8- 

{1/2Ay{Bid-'B^ - B,d-'Bl)qj>'j^q 

6 

V9+ 

{l/2A^)t''-'^%BldpB^ + B^.dpBDq'y^q 

6 

V9- 

{l/2A^)lt^''"‘%Bld,B,z - B.dpBDq'f^q 

6 

V10+ 

{l/2A^)e^'-'P^BtdpB,z + B,dpBl)q-fqj^q 

6 

VIO- 

{l/2A^)ie>'-'^‘^{BldpB,z - BodpBl)q-f^'y^q 

6 


TABLE I. Possible Hermitian contact operators up to dimension 6 that couple spin-1 dark matter to SM quarks (or other 
fermions). 


The operators in Table |T] are assumed to be the low energy manifestations of some more fundamental ultraviolet 
theory. All of the contact operators given above can arise from renormalizable interactions in which the dark matter 
pair is produced by the s-channel exchange of a spin-1 or spin-0 mediator (for example, a heavy Z' or Higgs particle), 
or by the A/n-channel exchange of a spin-1/2 particle. A detailed analysis of these constructions is presented in |18j . 
Our intention is to perform a completely general analysis of spin-1 dark matter at the LHC; under the assumption 
that the effective theory description is valid at these energies, the operators in Table |l] represent a basis set to describe 
these interactions. 

We are interested in the scaling of the matrix element for the process qq ^ B with respect to the energy E in 
the center-of-mass frame of the qq system. This scaling is determined by the following considerations: 

• The matrix element scales by a factor {E/A)'^~'^, where d is the dimension of the operator and A is the energy 
scale of the coefficient. 

• The matrix element also scales by additional factors of E/ms for each DM longitudinal polarization vector. 
The number of such vectors can be found by determining the C, P, and J quantum numbers of the DM state 
that can be created by each operator. These in turn determine the L and S quantum numbers of the DM final 
state, which determine the polarization vectors. 

The C and P quantum numbers of the dark matter two-particle state (for this purpose, the jet(s) in the final state 
are irrelevant) created by each operator are determined by the transformation properties of the DM bilinear. The J 
quantum number of the dark matter two-particle state is determined by the rotational transformation properties of 
the DM bilinear. The L and S quantum numbers of the DM final state system are then given by 

C:(-l)^+^, \L-S\<J<\L + S\, (1) 

which are valid when the dark matter is a boson. Following in Table |h] we display the quantum numbers of the 
DM state created by the various terms in the above operators. Note that we have ignored all operator terms involving 
the quark bilinears q^^q and qj^j^q. The former vanishes identically when acting on any quark-antiquark initial 
state, while the latter yields a matrix element that is proportional to rriq, and thus vanishes in the limit niq <C A, as 
is relevant here. 

Using Table [Til write the DM hnal state as a linear combination of states in the \L, Stot, J, Jz) basis 

(where Stot is the total spin of the DM system), then rewrite the state in the \L, Lz, Stot, Stot basis, and finally 
rewrite the state in the \L, Lz, Si, Siz, S 2 , S 2 z) basis (where Si and S 2 are the spins of each of the two dark matter 
particles.). The matrix element then receives a factor E/vtib enhancement for each DM particle with spin projection 
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Operator Term Cb Pb J State 


VI 


-b 

+ 

0 

L = 0, 5' = 0;L = 2, 5 = 2 

V2 

{l/A)zBlB'^q'f^q 

-b 

+ 

0 

L = 0, 5' = 0;L = 2, 5 = 2 

V3 

il/2AyiBtdiB‘' - B''diBl)q-i^q 

- 

- 

1 

L = 1, 5 = 0; L = 1,3, S' = 2 

V4 

{l/2A^)i{BldiB'' - B-'diBtW-z^q 

- 

- 

1 

L = 1, 5 = 0; L = 1,3, 5 = 2 

V5 

{l/A)iBlBjqa^^q 

- 

+ 

1 

L = 0,2, 5 = 1 


{1/2AXBIB. - BlBoq)qa°^q 

- 

- 

1 

L = 1, 5 = 0; L = 1,3, 5 = 2 

V6 

{l/A)BlB,qa^^y^q 

- 

+ 

1 

L = 0,2, 5 = 1 


(l/2A)(RtR-RjBo)g<TVg 

- 

- 

1 

L = 1, 5 = 0; L = 1,3, 5 = 2 

V7+ 

{l/2A^){Btd-'B, -b B,d''B])qYq 

-b 

- 

1 

L = 1, 5 = 1 

V7- 

{l/2A^XBld-'B, - B^d‘'Bl)qYq 

- 

- 

1 

L = 1, 5 = 0; L = 1,3, 5 = 2 

V8+ 

{l/2A^){Btd^B, -b B,d''Bl)qYy^q 

-b 

- 

1 

L = 1, 5 = 1 

V8- 

{l/2Ay{Btd‘'B, - B,d‘'Bl)qY'y^q 

- 

- 

1 

L = 1, 5 = 0; L = 1,3, 5 = 2 

V9+ 

{l/2A^)P°^'^{BldjBk + BodjBl)qjiq 

-b 

+ 

1 

L = 2, 5 = 2 

V9- 

{l/2A^)iP°^^{BldjBk - BodiBl)q~iiq 

- 

+ 

1 

L = 0,2, 5 = 1 


{l/2A^)iP^°\B]doBk - BjdoBl)q'^iq 

- 

+ 

1 

L = 0,2, 5 = 1 

V10+ 

{l/2A^)P°^\Bld,B^ + B,d,Bl)q^,^\ 

-b 

+ 

1 

L = 2, 5 = 2 

VIO- 

{l/2A^)P°i\BldjBu - BodjBl)q'yi-i\ 

- 

+ 

1 

L = 0,2, 5 = 1 


{l/2A^)P^°\BldoBk - B,doBl)q'ya^q 

- 

+ 

1 

L = 0,2, 5 = 1 


TABLE II. The charge conjugation (Cb), parity (Pb), and total angular momentum (J) quantum numbers of the DM system, 
as well as possible orbital (L) and spin (5) angular momenta of the DM state. Only nonvanishing or non-negligible terms in 
each operator are shown. 


’(L2)= 


= 0. The details of this derivation are provided in Appendix A, and we summarize the energy dependence 

While it may be possible that the coefficient of the leading term experiences an 


of the leading term in Table III 


Operators 

Dimension enhancement Polarization enhancement 

VI, V2, V5, V6 

E/A 

{E/ruBf 

V3, V4, V7_, V8_ 

{E/Af 

(E/mBf 

V7+, V8+, V9±, V10± 

{E/Af 

E/niB 


TABLE III. The energy enhancement factors in the matrix element for the process qq —>■ B. 

accidental cancellation, explicit calculation of the matrix elements indicates that this is not the case. The squared 
matrix elements are listed in Appendix B. 


III. CONSTRAINTS FROM UNITARITY 

If the scale of the new physics mediating the dark matter-quark interaction is sufficiently light, the contact- 
operator approximation will break down. In simple models with a single mediator, the energy suppression scale A of 
the contact operator is generally larger than the mass scale of the mediator, implying that one should not trust the 
contact-operator approximation for processes where the DM system center-of-mass energy is larger than A. But in 
more complicated models with large numbers of mediators, it is not clear that the mediator mass must be smaller 
than A. More generally, the tightest constraint one can rigorously impose is that the dark matter production matrix 
element satisfy unitarity when evaluated at the energy of the hard process. 

To impose this constraint, we follow the formalism and notation of m- We thus consider the matrix element for 
the on-shell process qq —>■ B. If the initial state is a helicity eigenstate, then the fundamental unitarity constraint 
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Operator 
VI, V2 
V3, V4 
V5, V6 
V7+,V8+ 
V9+,V10+ 
V7-, V8_ 
V9-, VIO- 


Constraint 


E-sjE'^ - m] 
E{E^ 


72-k‘2A* 


3 + ^{E^ - ml) ] <1 


3 + ^{E^ - ml) ] <1 


E^E^ - 


72.^A^ + 


£'3(^;2 _m|)3/2 
187r2mgA‘^ 

EjE'^ -ml)^^^ 

187r2mgA‘^ 


E'^E"^ - ml)'^/^ 
ISTT^m^A^ 

E^{E'^ - mlf/'^ 
327r2A'‘ 


1 + 


1 + 2 


< 1 

< 1 
E'^ 


E^ 


< 1 


< 1 


Benchmark Amin (TeV) 
1.59 X 10® 

274 

5.31 X 10'‘ 

8.66 

8.66 

274 

8.66 


TABLE IV. Unitarity constraints on the energy E of a dark matter particle in the center-of-mass frame of the qq —>■ B 
interaction. These constraints can be rephrased in terms of a maximum B^B invariant mass, which is then applied during 
event generation to get conservative collider sensitivities. We have also included the minimum value of A from these constraints 
for each operator, using the benchmark values of i? = 1 TeV and ms = 1 GeV. 


can be written as 




( 2 ) 


where we have expanded the matrix element M.i^f in Wigner d-functions as 

OO 

M.^f(B) = . (3) 

1=0 

Here, Xij are the total helicities of the initial and final states, respectively, j is the total angular momentum of the 
state, and 0 is the scattering angle. The Wigner d-functions which are relevant here are given by 


wo _ 1 

“0,0 ~ 3-, 


'^1,0 = “'^-i.o ~ '^0,-1 = ~^o.i = ~ 
<^1,-1 = = 


1 + cos 0 
2 

sind 
1 — cos 0 


dj 0 = cos 9. 


( 4 ) 


They satisfy 


dcosO dy^dy^ — ^ 

diy{e = 0) = Sxy. 


(5a) 

(5b) 


The unitarity constraints (on the qq ^ B process) shown in Table IV are obtained by explicitly evaluating Eq. © 
for each of our 14 operators. We also give the minimum value for A allowed by these constraints, using the benchmark 
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values of = 1 TeV and tob = 1 GeV. Our analysis will apply these constraints on an event-by-event basis: If we 
find that the LHC is sensitive to a particular value of A, the events used to establish this sensitivity had dark matter 
energies satisfying these constraints. Note that Eq. (H provides constraints for each initial quark-antiquark state for 
which the matrix element is nontrivial. For each of the operators there are at most two initial helicity eigenstates 
that are relevant. But for each operator with multiple possible initial helicity eigenstates, the constraints arising from 
each of those possible initial states are degenerate. 

In fact, because an energetic monojet is emitted, the actual hard process is 2 —> 3 with an off-shell intermediate 
(anti)quark. The unitarity bound above constrains the matrix element for the subprocess qq —>■ XX, where the quark 
and antiquark are on shell. However, for the energies and cuts relevant for the LHC analysis, the virtuality of the 
quark is a subleading effect, and a correct accounting for this virtuality results in only a small change in the unitarity 
constraint. One could also use the unitarity condition to directly constrain the matrix element for the full 2 —>■ 3 hard 
process, but this constraint is weaker than that arising from applying the unitarity condition to the 2 —> 2 subprocess. 
A more detailed discussion of these issues can be found in m- 


IV. LHC BOUNDS 


We consider a search for monojet signatures arising from the process pp —> B -1-jets, assuming that dark matter- 
quark interactions arise from any of the 14 contact operators shown in Table |Tj We assume a universal coupling to up 
and down quarks, and no coupling to the heavier generations. Signal and SM background events are generated using 
the MadGraph/Pythia/Delphes simulation chain. The 14 operators are input into MadGraphS, and simulated 
pp B^Bj events for each operator are generated using MadEvent [20j. These events are then showered and 
hadronized using Pythia-6.4 and event detection at ATLAS is simulated using Delphes-3 |22]. 

We impose the unitarity constraints in Table IV on an event-by-event basis [53]. For each generated event for 
which the qq ^ B^B matrix element satisfies Eq. (|2|), the contact-operator approximation can provide an adequate 
description of the physics. For events where Eq. ® is not satisfied, some new physics must come into play; to be 
conservative we simply reject those events. Note that the unitarity constraint for a specific choice of DM mass and A 
is entirely a function of the invariant mass of the final state dark matter system; since the final states of the selected 
events are different from the final states of the rejected events, interference effects are automatically removed. 

However, it is important to point out that for events that satisfy Eq. ([^, one only knows that the contact approx¬ 
imation is consistent with unitarity. Although it is not necessary for there to be any additional non-negligible new 
physics, any particular model may exhibit new physics not captured by the contact-operator approximation which is 
relevant at energies below the limit at which unitarity is violated. In fact, it might be surprising if new physics fixed 
unitarity for processes at energies where the contact approximation would slightly violate unitarity, but had only a 
negligible impact even at slightly lower energies where the contact approximation would satisfy unitarity. 

We also note that the unitarity analysis of m and used in this work assumes that the incoming states are in 
a helicity eigenstate. However, for actual events at the LHC, the incoming partons of the hard process will not 
generally be in a helicity eigenstate, but rather in a linear combination of helicity eigenstates. To generalize the 
unitarity analysis to a generic initial state is beyond the scope of this work, but we expect the maximum invariant 
mass for which unitarity is satisfied to change by at most an C>(-\/2) factor. Thus, although we have made largely 
conservative approximations, the above caveats suggest that our limits may be uncertain by factors ~ 0{\/2). 


A. Constraints from the 8 TeV LHC 

The ATLAS Collaboration has performed a search [24] for the monojet signal using 20.3 fb~^ of data at a center-of- 
mass energy of ^/s = 8 TeV. The main event selection criteria are the requirement of at least one jet with transverse 
momentum pT > 120 GeV, several jet isolation and quality-control requirements, and a large amount of missing 
energy. ATLAS chose several signal regions, defined by differing amounts of minimum missing energy, in order to 
tune their analysis to various sources of new physics. In this work, we will use their signal region SR4, requiring 
missing energy > 300 GeV. 

The ATLAS analysis specifies the full set of triggers, jet candidate requirements, preselection cuts, and final signal- 
region missing-energy cuts performed on the 20.3 fb“^ data set. For this analysis, we use a slightly simplified subset 
of these cuts which are appropriate for our Monte-Carlo-generated event simulation chain and which capture the 
primary features of the ATLAS selection. First, an initial cut of > 250 GeV is applied at the parton level in order 
to decrease the event veto rate and increase statistics after final cuts. The detector level events are then required 
to satisfy the selection criteria given in Table [Vj We have tested these cuts by applying them to Z —>■ vv, W vl. 
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Jet reconstruction: 

anti-fer, using R — 0.4 

Jet definition: 

Pt > 30 GeV and It;! < 4.5 

Lepton veto: 

electrons: pr > 7 GeV and I 77 I < 2.47 
muons: pr > 7 GeV and \rj\ < 2.5 

Leading jet: 

PTji > 120 GeV and |» 7 | < 2 
PTjilt^T > 0.5 

Separation (all jets): 

Acj>{pTj,$T) >1-0 

Missing energy: 

$T > 300 GeV 


TABLE V. Monojet selection cuts for the ^/s = 8 TeV LHC analysis. Acj) is the angular separation between the selected jet 
and the missing transverse momentum, R is the radius parameter used in the anti-fer jet clustering algorithm |25| . and r; is the 
pseudorapidity. 


10® 


% 10 « 
< 

10 ^ 


102 

1 10 100 500 



uiB (GeV) 

FIG. 1. The ATLAS 95% C.L. exclusion bounds on the vector DM contact-operator scale A, using the 20.3 fb“^ data set 
at -^i = 8 TeV. All even-numbered operators (V2, V4, ...) are visually indistinguishable from their odd-numbered operator 
counterparts, and thus are not. The shaded region at A < 2ms represents the regime where the effective-operator description is 
naively expected to break down. Note that all events used to establish these sensitivities were required to satisfy the unitarity 
constraints on four-point interaction energies. 


and W ^ i^T background events produced using our same MadGraph/Pythia/Delphes simulation chain, and find 
that our predicted event rates match the ATLAS event rates to within ~ 5%. 

Using the SR4 signal region, ATLAS is able to exclude at the 95% C.L. any new-physics source of monojet events 
which gives rise to a cross section of 51 fb or greater. We use this constraint to bound the new physics scales A for 
each of the 14 vector DM contact operators. We apply the set of kinematical cuts in Table [V| to hnd the total event 
rate (or, the total cross section after cuts ctdm) for a given A. 

To impose unitarity constraints, as discussed above, we apply at an event-by-event level a cut on the maximum 
DM invariant mass, or equivalently, on the center-of-mass energy of the underlying DM-SM four-point interaction 
^sdm- Applying this cut safely underestimates the total event rate while excluding events from regions of phase space 
where unknown high-energy physics is required by unitarity. With this invariant mass cut in place, we tune A so that 
the total event rate after all cuts corresponds to the new-physics cross section excluded by ATLAS, ctdm = 51 fb. 
This provides us with 95% C.L. exclusion bounds on A for each operator, and over a range of DM masses ms; these 
exclusion bounds are shown in Fig. We note here that the exclusion bounds for the even-numbered operators (V2, 
V4, ...) are visually indistinguishable from their odd-numbered counterparts over the mass range of Fig.[^ and hence 
are not included in the figure. 
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Jet reconstruction: 

anti-fer, using R = 0.4 

Jet definition: 

Pt > 50 GeV and \r]\ < 3.6 

Lepton veto: 

electrons: pr > 7 GeV and \ri\ < 2.47 
muons: pr > 7 GeV and \rj\ < 2.5 

Leading jet: 

PTji > 300 GeV and \p\ < 2.0 
PTjilt^T > 0.5 

Separation (all jets): 

l\cj>{pTj,$T) > 1-0 

Missing energy: 

$T > 600 GeV 


TABLE VI. Monojet event selection cuts for the ^/s = 14 TeV LHC analysis. These cuts are similar to those used for the 
8 TeV LHC, with adjustments to jet definition, leading-jet momentum, and missing energy, motivated by the higher energy 
and pileup at 14 TeV |27) . 


Operators V3, V5, V 6 , V7, and V9 would also permit the decay of a = 1 quarkonium state to B, if 
kinematically allowed. Constraints on these operators from bounds on the invisible decay rate of T(IS') were found 
in [26) . in the case where dark matter couples to 6 -quarks. We can compare those results (assuming a coupling to 
6 -quarks) to the LHC sensitivity found above (assuming instead a coupling to first-generation quarks). In the region 
of overlap, the LHC reach in A is at least a few orders of magnitude larger than that obtained from current bounds 
on invisible quarkonium decay. 


B. Sensitivity of the 14 TeV LHC 


With the recent commencement of the high-energy LHC run, a new data set at increased center-of-mass energy 
will soon be in hand. We now calculate the expected sensitivity of this new data set to spin-1 dark matter, coupling 
to Standard Model quarks through these 14 operators. In order to best compare with the current exclusions found 
above, we calculate the expected 95% C.L. exclusion bounds under the assumption that future data sets do not see 
any events above the expected background. For this analysis, we will assume that the LHC Collaboration will soon 
increase beam energies to their full design energies for a total center-of-mass energy of ^/s = 14 TeV, though these 
results should still be insightful if the LHC remains operating at yfs = 13 TeV. 

To calculate exclusion bounds for the 14 TeV run, we must first know the event rate of SM irreducible background 
events at this energy. We estimate this event rate by simulating the processes Z —>• w, W —and W —>■ lyr 
using the same Monte Carlo simulation chain and selection cuts used to calculate the signal event rate. According to 
previous monojet analyses, including the ATLAS 20.3 fb“^ analysis [24], these three processes contribute ^ 95% of the 
total irreducible background to the monojet signal, so we neglect other sources such as tt and single top production. 

The simulation chain is identical to that used for the 8 TeV data set described above with adjustments to the 
parton-level and detector-level jet dehnitions and cuts. These adjustments are motivated by an ATLAS Collaboration 
study of monojet searches at 14 TeV m- Jets are defined to have a higher minimum transverse momentum of 
Pt > 50 GeV in order to mitigate the higher pileup expected at this increased energy. Additionally, all jets are 
required to be more central, with I 77 I < 3.6, as jet reconstruction algorithm performance is not well understood in 
the forward regions. The leading-jet minimum px and minimum missing energy are both increased to reflect the 
increase in available energy to pxji > 300 GeV and > 600 GeV, and the parton-level cut on missing e nerg y is 


VI 


correspondingly increased to > 500 GeV to increase statistics. All hnal selection cuts are shown in Table 

On applying these cuts to the generated background events, we find the cross sections for the three dominant 
background processes, as shown in Table |VlT| These event rates are represented by usm x e, where ctsm represents the 
cross section for the process as calculated from the MadGraph/Pythia/Delphes simulation chain for our choice of 
parton-level cuts, and e represents the efficiency of the selection cuts given in Table jVTj The product ctsm x e is the 
physically meaningful quantity, as both ctsm and e depend on our choice of parton-level cuts. 

Using these results, we calculate the total new-physics cross section that can be excluded by the 14 TeV LHG for 
a given integrated luminosity. The 95% G.L. expected exclusion bound is given by 




A’dm(A)^ 


-^SM ' 


■ A^dm(A) -I- Z? 


= 3.84, 


( 6 ) 


SM 


where A^sm and AAim are the number of expected background and signal events, respectively, and Ssm is the systematic 
uncertainty in the number of background events. We can rewrite the above condition in terms of integrated luminosity 










CTSM X e (fb) 

Z ^ vv 

79.93 

W —>■ TVr 

16.99 

W Ivi 

10.51 

Total 

107.4 


TABLE VII. Background cross sections at the ^/s = 14 TeV LHC for the three dominant background processes using the 
cuts given in Table |Vl| As determined by previous monojet analyses, these three processes contribute ~ 95% of the total SM 
irreducible background, so we neglect other sources in our background estimation. 



niB (GeV) 

FIG. 2. The expected ATLAS 95% C.L. exclusion bounds on the vector DM contact-operator scale A, using 100 fb~^ of data 
at y's = 14 TeV. All even-numbered operators (V2, V4, ...) are visually indistinguishable from the odd-numbered operators 
shown. The shaded region at A < 2m b represents the regime where the effective operator description is naively expected to 
break down. 


as 


ctdm(A)^ • C 


CSM + crDM(A) 


■ ^Im 


'SM 


c 


= 3.84, 


(7) 


where cr = tr x e is the cross section after cuts, and Jsm is the percent systematic uncertainty in the number of 
background events. As determined by the ATLAS Collaboration, the uncertainty in background events at 8 TeV 
ranges from ^ 2% to 10% (see, for instance. Tables 4 and 5 of [21]), so for this analysis we take (5 sm = 0.05. We note 
here that although the LHC sensitivity to each individual operator will depend on this uncertainty, our plots which 
overlay all vector operators extend over such a large range that small variations of (5 sm will appear negligible. For this 
analysis we take C — 100 fb“^, which is the approximate integrated luminosity expected to be collected by the end o f 
the first LHC run at 14 TeV. For these values of £, i5sm, and (tsm (as determined above and displayed in Table VH|, 


we find that the LHC at 14 TeV is able to exclude at the 95% C.L. a new-physics event rate of crDM(A) = 10.74 fb. 
This sensitivity is dominated by systematic uncertainties, and thus roughly scales linearly with 5 sm- 

Analysis of the signal events for each of the 14 contact operators proceeds in the same way as before, where A is 
tuned such that the total event rate after cuts, including the cut on DM invariant mass constrained by unitarity, is 
such that crDM(A) = 10.74 fb. The exclusion bounds on A over a range of DM masses are shown in Fig. Once 
again, only the odd-numbered operators are shown, as the even-numbered operators are visually indistinguishable 
from their odd-numbered counterparts on this plot. 
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A oc 


VI, V2, V5, V6 


V3, V4, V7-, V8- 


V7+, V8+, V9±, V10± 




A4m| 


E 




E 

ruB 

~E 

ms 

~E 


TABLE VIII. Scaling of the qq —>■ B cross section and overall collider sensitivity to A based on the enhancements from 
operator dimension and longitudinal polarizations given in Table [111] We see that the three subsets of operators should scale 
differently, which is indeed the behavior seen in Figs. and 


C. Discussion of sensitivities 


The results shown in Figs. and [^encode the enhancements each operator receives from both operator dimension 
and longitudinal polarizations, and the features in these plots can be fully understood by a simple accounting of 
these enhancements. From Table m we see that there are three subsets of operators that are defined by how 
they scale: operators V(l,2,5,6) scale as E^/{Am‘^), operators V(3,4,7_,8_) scale as and operators 

V(7+,8+ ,9±, 10±) scale as E^/{A?mB)- Thus depending on the operator there will be three different ways in which 
the pp —>■ B cross section, and thus the collider sensitivity to A, scales. This is shown in Table VIII and the 
behavior can clearly be seen in the LHC exclusion bounds shown in Figs. [T] an d [2} 

Furthermore, the overall magnitude of the exclusion limits in Figs, [^and^ matches what we expect. For instance, 
at the LHC we expect the center-of-mass energy of these quark-D M inter actions to be 0(1 TeV); given a DM mass 
of tob = 1 GeV, we then expect sjE/nriB 


VHI 


the sensitivities to A for the three different 


30. According to Table 

subsets of operators should then, for ms = 1 GeV, stand approximately in the ratio 1 : 30 : 30000. This is indeed 
what is seen in the figures. Thus it becomes clear that collider sensitivity to vector-DM production can receive an 
enhancement of several orders of magnitude at the LHC simply due to the presence of the longitudinal polarization 
mode. We now study how these collider bounds on A translate to constraints on scattering cross sections at direct- 
detection experiments. 


D. LHC bounds on velocity-independent scattering 


Operators VI and V3 permit velocity-independent spin-independent (SI) scattering, while operators V5 and VI0+ 
permit velocity-independent spin-dependent (SD) scattering (see, for example, [H]). For these operators, LHC con¬ 
straints on the energy scale of the operator (A) can be expressed as constraints on the DM-nucleon scattering cross 
section. 

The DM-nucleon scattering cross sections may be written as 




IJ-N 


\Q-nm?grnj^i 


E i-^i' 


( 8 ) 


where /iy is the reduced mass of the DM-nucleon system, rrijv is the nucleon mass, and M. is the scattering matrix 
element. Assuming isospin-invariant couplings to first-generation quarks, the scattering cross sections can be written 
as (see also [2H] 1 


Af(Vl) 

^SI 


N(V3) 

o-si 


Af(V5) 

^SD 


iV(V10+) 

^SD 




Lttto^A^ 




l^N 

dTrA"^ 


2'!Trni,A'^ 


27rA4 


(^B^ipv) ^ 


N(v) 


y 


N{t) 


B 


,N{pv)' 

d 




(9a) 

(9b) 

(9c) 

(9d) 
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FIG. 3. The ATLAS 95% C.L. exclusion bounds on using the 20.3 fb“^ data set at y/s = 8 TeV, as well as 

the expected 95% C.L. sensitivity curve for a 100 fb“^ data set at y/s = 14 TeV. Also plotted are exclusion bounds from 
SuperCDMS [32], CDMSlite [33] . and LUX (31], as well as a representative minimum cross section below which neutrino- 
nucleus coherent scattering is significant |36 |. 


where are the nucleon form factors for u- and d-quarks for scalar, vector, pseudovector, and tensor 

structures. The vector nucleon form factors are fixed by gauge invariance, and are given by 


DP(v) _ d"(*') 


= 2 , 


B 


n(v) 


= B 


pM 


= 1 . 


( 10 ) 


The remaining nucleon form factors are subject to uncertainties related to the structure of the nucleon. We use the 
following values as benchmarks [^51t30j : 


RP(s) , 

ryn(s) 

^10, 

R"(«) , 

nPi^) 

-7, 


RP(PJ') ^ 

Tyn{pv) 

^d 

- 0.84, 

Tjn(pv) ^ 

dP(p«) 

^d 

- -0.43, 


fjpW , 

rjn{t) 

- 0.54, 


nPit) 
Bd - 

- -0.23. 

(11) 


We plot current LHC 95% C.L. bounds on crsi and ctsd in Figs. and respectively. We also plot current 
direct-detection limits from LUX |3T|, SuperCDMS [35], and CDMSlite [33] for spin-independent scattering, and 
PICO-2L [33] and XENONIOO [35] for spin-dependent scattering. In addition, we plot the minimum spin-independent 
scattering cross section that can be probed while still having an insignificant contribution of events from neutrino- 
nucleus coherent scattering [36] . 

For operators that produce spin-independent scattering, the LHC is only competitive with direct-detection exper¬ 
iments at low mass. Note that for operator VI, although the LHC energy reach A increases dramatically for small 
tub, its sensitivity to trsi does not increase as dramatically because oc Ijm^. Additionally, despite the fact 

that collider sensitivity to A for VI far exceeds the sensitivity for V3, collider bounds more tightly constrain ctsi for 
V3. 

For spin-dependent scattering, however, the LHC sensitivity far exceeds that of current direct-detection experiments. 
Additionally, the LHC is more sensitive to vector operators V5 and V10+ than it is to fermionic operators which lead 
to SD scattering, such as xi^l^X and This is also the case for vector and fermionic spin- 

independent operators. 
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tub (GeV) 

FIG. 4. The ATLAS 95% C.L. exclusion bounds on using the 20.3 fb“^ data set at y/s = 8 TeV. Also plotted are 

95% C.L. exclusion bounds from P1CO-2L |34| and XENONIOO |35| . 


V. CONCLUSION 


We considered the sensitivity of the LHC to spin-1 dark matter that couples to first-generation quarks via a contact 
operator. We found that a monojet search strategy can probe contact operators with an energy suppression scale A 
which can range from 0(1 — 10®) TeV. This large energy reach results from dark matter final states which include 
one or more longitudinally polarized spin-1 particles, yielding large enhancements to the production matrix element. 
In this analysis, we imposed conservative cuts to ensure that we only consider regions of phase space for which 
the contact-operator approximation is consistent with unitarity. For operators that permit velocity-independent SD 
scattering, these bounds far exceed those obtainable from current direct-detection experiments, while for operators 
that permit velocity-independent SI scattering, these bounds are comparable to those obtainable from direct-detection 
experiments for light dark matter. 

It is worth noting the implications of these results for future high-energy hadron colliders. As indicated in Table [VlITl 
the sensitivity of monojet searches to spin-1 dark matter scales as a high positive power of the characteristic energy 
of the collider; for operators V(l, 2, 5, 6), the sensitivity scales as A oc £1® and an increase in collider energy is much 
more effective than an increase in luminosity. Indeed, for these operators, the sensitivity of a 14 TeV run of the LHC 
is approximately an order of magnitude greater than that of an 8 TeV run. One might therefore expect that a future 
0(100) TeV hadron collider could provide a sensitivity orders of magnitude greater than current bounds. 
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Appendix A: Energy enhancement due to longitudinal polarization 


The two-particle spin state in the [S'!, S 2 , Stot, Stot^) basis can be expressed in the 151, Siz] 82 , S 2 z) basis as follows: 


| 1 , 1 , 2 , 2 ) 
| 1 , 1 , 2 , 1 ) 

|1,1,2,0) 

|1,1,2,-1) 

|l,l,2,-2) 

|1,1,1,1) 

|1,1,1,0) 

|1,1,1,-1) 

|1,1,0,0) 


|i,i;i,i), 

^|1,0;1,1) + ^|1,1;1,0), 

^| 1 ,- 1 ; 1 , 1 ) + ^| 1 , 0 ; 1 , 0 )+ ^| 1 , 1 ; 1 ,- 1 ), 
^| 1 , 0 ; 1 ,- 1 ) + ^| 1 ,- 1 ; 1 , 0 ), 


|i,-i;i,-i), 

^| 1 , 0 ; 1 , 1 )-^| 1 , 1 ; 1 , 0 ), 
^|i,-i;i,i)-^|i,i;i,-i), 
^| 1 , 0 ; 1 ,- 1 )-^| 1 ,- 1 ; 1 , 0 ), 
^|i,-i;i,i)-^|i,0;i,o) + ^|i,i;i,-i). 


(Al) 


Since the matrix element receives an A/tob enhancement factor for each longitudinal polarization, we find the following 
relations between the DM total spin state and the leading energy enhancement: 

{E/msf : Stot = 2, Stot^ = 0; Stot = 0, Stot^ = 0 , 

{E/mBf : Stot = 2 , Stou = ± 1 ; Stot = 1 , Stou = ± 1 , 

{E/mBf : Stot = 2, Stou = ±2; Stot = 1, Stot^ = 0 . (A2) 


The last step is simply to consider each |L, Stot, J, Jz) final state allowed by every term of every operator [T7], and 
expand in the basis |L, Lz,Stot, StotJ via Clebsch-Gordan coefficients. We then find the leading energy enhancement 
allowed, yielding 

L = 0,Stot = 0,J = 0^ (E/mB? , 

L = 2, = 2, J = 0 ^ (EjmBf , 

L = 1, = 0, J = 1 ^ (A/mB)^ 

L = 0, Stot = 1,,/ = 1 —>■ E/niB , 

L = 1, Stot = 1, T = 1 —> E/iTiB , 

L = 2, Stot = 1, T = 1 —> E/mB , 

L = 2, = 2, J = 1 ^ E/ruB ■ (A3) 


Note that these expressions are sufficient to compute the leading energy enhancement for each of the operators we 
consider. 


Appendix B: Squared matrix elements 


Here we list the squared matrix elements (summed over spins and polarizations) for the process qq —>■ B^B. 
We take E to be the energy of a DM particle in the center-of-mass frame of the dark matter system, so that 
E'^ = k'^ +iTiq = kg The angle between the axis of the quark-antiquark system and the axis of the dark matter 

system is 9. In the limit > 0, we find 


|AfviP= X! l-^V2p 

spins spins 

|Afv3p = l-^V4p 

spins spins 



SE'^Pg 

A4 


sin^ 9 




(Bl) 

(B2) 



13 


iAiv6p 

spins spins 

iaiv7+p= y^ iaiv8+p 

spins spins 

y^ l■A^V9+P= X! I-^VIO+P 

spins spins 


IG-E^ / f fcn 2 /I ' 

—^ l + 2-|- 1 + -f-cos^e 

\ \ mi , 


A4 ^2 (1 


COS' 


m 


s 




16E^k% 
m\ ^ 


COS' 


^0). 


For the V(7 — 10)_ operators, we find 

y^ |Adv7_|^= X] |A1 v8_|^ 

spins spins 

y^ iAdv9„p= y]] iAivio_p 

spins spins 


16E^ fc| 

16£;4 


1 + cos^ 9 + 2—^ sin^ 6 

mi. 


(l + cos^ O) + sin^ 9^ . 


(B3) 

(B4) 

(B5) 


(B6) 

(B7) 


Appendix C: Properties of operators V(7— 10)_ 

Operators V(7 — 10)_ were not discussed in [T7]. We repeat much of that analysis for these operators, for com¬ 
pleteness. In Table [P^ we write the annihilation matrix element factors that arise from the DM bilinears relevant for 
operators V(7 — 1^_. We assume that the DM system is in the center-of-mass frame, with B {B^) having spatial 
momentum k {— k). In Table |xj we write the scattering matrix element factors in the center-of-mass frame that 
arise from the same bilinears, where the incoming DM particle has polarization vector e and the outgoing particle 
has polarization vector e'. In Table pOl we indicate if the operator permits DM s-wave annihilation, and indicate 
the factors of momentum transfer (q) and of DM velocity perpendicular to momentum transfer (u'*") that suppress 
the spin-independent and spin-dependent scattering cross sections (including if the factor arises from the DM or SM 
bilinear). 


Bilinear 

Annihilation matrix element 

- B.,d,Bl) 

- BjdoBl) 

-eoijk{BldjBk - BodjBl) 
-ekijo{BldjBo - BkdjBl) 
B^''d,Bo - B’'d,Bl 
B^''d,B, - B''d,Bl 

0 

itijkE 

—leijhE + £ 2 ^ 1 ) 

0 

2%E{t%e\ — 6162 ) 


TABLE IX. The annihilation matrix element factors for spin-1 dark matter bilinears relevant for operators V(7 — 10) _. 


Bilinear 

Scattering matrix element (SD) 

(Bld’^B^ - B.d-^Bl) 


e°''‘^'^{BldpB^-B^dpBl) 


e*''‘=^BldpB, - B^dpBl) 

0 


TABLE X. The scattering matrix element factors for spin-1 dark matter bilinears relevant for operators V(7 — 10) _. 
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Operator 

Structure 

(Tsi suppression 

(TSD suppression 

s-wave? 

V7_ 


(SM); (DM) 

(DM) 

9^ (SM); q^ (DM) 

No 

V8_ 

(Bld-'B,, - B„d-'Bl)qj'^'f^q 

q'^v^^ (SM); q^ (DM) 

q^ (DM) 

No 

V9_ 

^e^'P<^(^BldpB„ - BpdpBl)qjp,q 

q^ (DM) 

q'^v^^ (SM); (DM) 

Yes 

V10 

^ppp-i^BldpBp - BpdpBl)q'yp^’^q 

0 

(SM); q^ (DM) 

Yes 


TABLE XI. The momentum (g) or velocity (w^) suppression of spin-independent or spin-dependent scattering cross sections 
mediated by operators V(7 — 10)Each suppression is labeled based on whether it arises from the Standard Model (SM) or 
dark matter (DM) bilinear. Also indicated is if the operator permits s-wave annihilation. If a cross section contains multiple 
terms with different kinematic suppressions, then they are listed on separate lines. 
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